Vibrational modes and lattice distortion of a nitrogen-vacancy center in diamond from 
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We investigate vibrational properties and lattice distortion of negatively charged nitrogen- vacancy 
(NV) center in diamond. Using the first-principles electronic structure calculations, we show that 
the presence of NV center leads to appearance of a large number of quasilocalized vibrational modes 
(qLVMs) with different degree of localization. The vibration patterns and the symmetries of the 
qLVMs are presented and analyzed in detail for both ground and excited orbital states of the 
NV center. We find that in the high-symmetry (C3„) excited orbital state a pair of degenerate 
qLVMs becomes unstable, and the stable excited state has lower (Cih) symmetry. This is a direct 
indication of the Jahn- Teller effect, and our studies suggest that dynamical Jahn- Teller effect in the 
weak coupling regime takes place. We have also performed a detailed comparison of our results with 
the available experimental data on the vibrations involved in optical emission/absorption of the 
NV centers. We have directly demonstrated that, among other modes, the qLVMs crucially impact 
the optical properties of the NV centers in diamond, and identified the most important groups of 
qLVMs. Our results are important for deeper understanding of the optical properties and the orbital 
relaxation associated with lattice vibrations of the NV centers. 

PACS numbers; 71. 55.-1, 71.15.Mb, 63.20.Pw, 78.55.-m 
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I. INTRODUCTION 

The negatively charged nitrogen- vacancy (NV) impu- 
rity centers in diamond have attracted much interest in 
recent years. They demonstrate a uniquely favorable 
combination of spin and optical properties: a NV elec- 
tronic spin has long coherence timei"— , the spin state 
of a single NV center can be optically initialized and 
read out^i^, and can be manipulated both optically2r— 
and magneticallyiSiiii. As a result, NV centers constitute 
promising candidates for applications in quantum infor- 
mation processing^^^^i, high-sensitivity magnetometry 
with nanoscale resolutioni^"— , and photonicsi^Ti^^. The 
NV centers also present an excellent platform for study- 
ing such fundamental problems of quantum mechan- 
ics as the dynamics of quantum spins coupled to their 
environmentii^i^^— , and for exploring quantum control 
and dynamical decoupling of solid-state spinsiii^"— . 

The electronic, optical, and spin properties of the NV 
centers are strongly affected by the lattice vibrations. 
For instance, due to slightly different atomic arrange- 
ment in the ground state and the excited state of the 
NV center, the optical emission and absorption is ac- 
companied by the lattice motion, i.e. involves the vibra- 
tional lattice excitation32r— . The impact of lattice mo- 
tion on the properties of the NV centers has been un- 
der investigation since 1970s22i^ till now^^i^^i"— . The 
lattice vibrations are responsible for the broadening of 
the light emission spectrum and appearance of additional 
spectral feature o^°'^^ , for depolarization of the emitted 
photons2ii^i^, and for the orbital and spin relaxation of 
the NV center o^^i'^^" — . The localized and quasi-localized 
vibrational modes are especially important for defect cen- 
ters in semiconductors, and are often used for investiga- 



tion of their properties'^""—. The detailed studies of the 
lattice vibrations, and of the (quasi) localized modes in 
particular, would be very useful for deeper understand- 
ing of the optical properties of NV centers, and would 
greatly improve our knowledge of the phonon-assisted 
spin and orbital relaxation of this system. Such stud- 
ies are also useful for improvement of the optical prop- 
erties of NV centers, as required for photonic-related 
applications^ii^""— . Very recently, an important first 
study of several localized modes has been reported for 
the NV center—. But detailed studies of other numer- 
ous localized and quasi-localized modes, which could be 
important for better understanding of the experimental 
data'^*'-'^^, are still lacking. 

In this paper, we present a detailed first-principle study 
of the localized distortions and vibrational properties 
of the negatively charge NV centers in diamond in its 
ground and excited orbital states. We show that the 
presence of the NV impurity noticeably modifies the vi- 
brational spectrum, and leads to appearance of numerous 
quasi-localized vibrational modes (qLVM) with different 
degree of localization. We present detailed description 
of several important qLVM, and analyze the changes in 
these modes upon optical excitation. Comparing our re- 
sults with earlier experiment a^"'^^ , we show that qLVM 
crucially influence the luminescence properties of the NV 
centers. The partial density of states (DOS) of the vi- 
brations (quasi)localized on the NV center is in excellent 
quantitative agreement with the results of Ref. [s^. This 
is the first, to our knowledge, first-principles explanation 
of these experiments. Our results also suggest that the 
one-phonon sidebands in the NV center emission and ab- 
sorption spectra are controlled by qLVM. We give a de- 
tailed description of the relevant modes, and show that 
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their properties (symmetry and energy) are in agreement 
with the experimental results^. 

Another interesting finding reported here is that some 
of the locahzed vibrational modes are unstable for the 
high-symmetry (Cs^) excited orbital state of the NV cen- 
ter. This is a direct indication that the Jahn- Teller (JT) 
effect influences the excited state, either by lowering the 
symmetry of the system (static JT effect), or by entan- 
gling vibrational and electronic states (dynamic JT ef- 
fect). Our calculations, as well as the known experimen- 
tal facts, evidence the latter. According to our results, 
the frequencies of the relevant modes are much larger 
than the JT energy, indicating the regime of the weak vi- 
bronic coupling, in agreement with known facts (such as 
e.g. small renormalization of the (/-factor and the shape 
of the optical transition bands). However, detailed quan- 
titative theory of the JT effect in NV centers requires a 
separate focused investigation, which is to be performed 
in the future. 

The rest of the paper is organized as follows. In Sec. II 
we describe the computational method used in the paper. 
In Sec. Ill we give a detailed description of the electronic 
structure of the negatively charged NV color center, and 
optimization of the atomic geometry. We discuss the lo- 
cal distortion of the atomic configuration in the ground 
and the excited states, and discuss the possibility of the 
static and dynamic JT effects in the excited state. In 
Sec. IV we present the results for the lattice vibrations 
in the ground and in the excited state. We discuss the 
phonon DOS, localization of different modes, and their 
properties. In Sec. V, we perform and extensive com- 
parison of our results with the known experimental data. 
Conclusions are given in Sec VI. 



II. COMPUTATIONAL METHOD 

The first-principles calculations were performed within 
the density functional theory (DFT) under the general- 
ized gradient approximation (GGA) of Perdew, Burke, 
and Ernzerhof (PBE)*'^, including spin polarization. The 
Vienna Ab-initio Simulation Package (VASP)^^ has been 
used to perform geometry optimization and force calcu- 
lations. A single electron was added to the electron oc- 
cupancy for the negatively charged vacancy. The extra 
electron is balanced by a uniform positive background to 
ensure the overall supercell is electrically neutral. We 
used F point for /c-point sampling and a plane wave basis 
set with a cutoff of 420 eV. 

The atomic structures of the NV center in its electronic 
ground and excited states have been optimized before 
the calculations for the vibration modes are performed. 
In the geometry optimization calculations, a 215-atom 
(3x3x3) cubic supercell with the box length of 3ao = 
10.719 A was used, where oq is the lattice constant of the 
diamond structure optimized by the VASP calculation. 
All internal atomic positions were relaxed until the forces 
were smaller than lO^'^ eV/A. 



The vibrational modes for the ground state and ex- 
cited state of the NV center were calculated within the 
harmonic approximation by first-principles calculations. 
The force-constant matrix of the supercell containing 
a negatively charged NV center defect are evaluated 
through the forces calculated via the Hellmann-Feynman 
theorem when each atom in the supercell is displaced one- 
by-one from its equilibrium position along ±x, iy, and 
±z directions respectively with a small displacement^^. 
We chose a displacement amplitude of 0.02 A from the 
atomic equilibrium positions in three Cartesian directions 
in the force calculation. The amplitude (0.02 A) is chosen 
here because it is small enough to remain in the harmonic 
part of the potential but large enough to avoid numer- 
ical noise. The forces acting on all other atoms within 
the supercell due to the one-by-one displacements of the 
atoms are used to construct the force constant matrix of 
the system. The 642 phonon frequencies plus three trans- 
lational modes and their corresponding eigenvectors are 
obtained by diagonalizing the 645 x 645 force constant 
matrix. 



III. ELECTRONIC STRUCTURE OF THE NV 
CENTER AND LOCAL DISTORTIONS FROM 
FIRST-PRINCIPLE GEOMETRY 
OPTIMIZATION 

The nitrogen- vacancy (NV) color center consists of a 
substitutional nitrogen atom (N) and an adjacent va- 
cancy (V) in the diamond lattice. The main features of 
its atomic and electronic structure are known rather well, 
and have been recently reviewed in detail in Ref.l43l For a 
negatively charged NV center, the six excess electrons are 
localized around the center. In the orbital ground state, 
the atomic arrangement has C'^y point symmetry, where 
one of the (111) crystallographic axes (passing through 
the N atom and the vacancy site) is the symmetry axis. 
Our calculations of the relaxed atomic structure are in 
agreement with these facts, see Fig. [TJ The N and C 
atoms at the nearest-neighbor shell of the vacancy re- 
lax outwards from vacancy. The atomic displacements 
from initial unrelaxed structure for the N and C atoms 
are 0.148 A and 0.101 A, respectively, and the distances 
from the vacancy to nitrogen and carbon atoms in the 
relaxed structure are 1.695 A and 1.646 A, respectively. 
These result s ag ree with those obtained previously, see 
e.g. Refs. H-lil 

The C^y point group has two types of non-degenerate 
singlet representations {Ai and A2), and a doubly de- 
generate representation E; below we will use lower-case 
letters to denote representations for single-electron states 
and vibration modes, and upper-case for many-electron 
states. The many-electron ground state of a NV center 
is an orbital singlet with the total electron spin 5* = 1, 
corresponding to representation (spin triplet, orbital 
singlet of A2 type). This state can be visualized within 
the qualitative model^'^ where six electrons localized at 
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FIG. 1: (Color online). Atomic displacements after geometry 
optimization, side view (top panels) and top view (bottom 
panels) with respect to the (111) axis. The red arrows show 
the atomic displacements, the arrow length is proportional to 
the diplacement magnitude. (a),(b): ground state, Csv sym- 
metry; (c),(d): excited state with Csv symmetry; (e),(f): 
the same excited state, but fully relaxed, with Cih symme- 
try. Yellow spheres denote N atom, blue spheres denote the 
C atoms adjacent to the vacancy, green spheres denote the C 
atoms adjacent to the N atom. 



FIG. 2: (Color online). The single-electron orbitals in the 
ground '^A2 and excited orbital states of the negatively 
charged NV center. The energies are given with respect to 
the top of the valence band. The states within the valence 
band are not shown. The symbols with the bar denote the 
spin-down states, the symbols without the bar correspond to 
the spin-up states. The states f and have ai symmetry, all 
other states have e symmetry. We use the notations of the 
parent Csv group to classify the states of the excited state 
with Cih symmetry, in spite of the reduced symmetry of the 
latter. 



the NV center populate the single-electron spin-resolved 
orbitals^i^, as shown in Fig. [5J Two electrons occupy 
the spin-up and spin-down levels immersed in the valence 
band, we exclude them from consideration. The next two 
electrons occupy the states v and D (spin-up and spin- 
down, respectively) of ai symmetry. These orbitals are 
in the band gap of the bulk diamond, and are localized at 
the carbons surrounding the vacancy and near the nitro- 
gen. The two remaining electrons occupy two degenerate 
spin-up orbitals of e symmetry {cx and Cy); the corre- 
sponding spin-down orbitals are unoccupied. The ener- 
gies of the single-electron states calculated here are in 
agreement with the previously published results^"— i^. 

Upon optical excitation the NV center undergoes tran- 
sition to an excited state of '^E symmetry (spin triplet, 
orbital doublet), which can be visualized as a transition 
of an electron from D to one of the spin-down orbitals 
ex or By. Starting from the optimized C^y geometry of 
the ground state, we found that the system conserves 
this symmetry after geometry optimization in the excited 
state. The outward atomic displacements of the N and 
C atoms are 0.090 A and 0.156 A, respectively, see Fig.[T] 
(middle panel), and the resulting C-to- vacancy and N- 
to-vacancy distances are 1.68 A and 1.63 A, respectively, 
also consistent with previous work^. 

However, if the geometry optimization for the excited 
state starts from a low-symmetry atomic arrangement, 
then the fully relaxed structure have the symmetry of 
the point group Cih- The signs of an instability of the 
Csy configuration are seen in the phonon spectrum: the 
high-symmetry Csy excited state possesses two degen- 



erate soft unstable phonon modes with formally nega- 
tive frequencies. The energy of the Cih atomic config- 
uration is 6.4 meV lower than the energy of the high- 
symmetry C31, configuration. While the displacements of 
the N atom in the C'lh and C^y configurations are the 
same, the displacements of the carbons adjacent to the 
vacancy are different: one moves by 0.147 A away from 
the vacancy site, and the two others move by 0.158 A, 
see Fig. [1] To check the possible finite-size effect in the 
geometry relaxation, we repeated this calculation for a 
twice larger supercell, with 511 atoms. The same result 
has been obtained, with the fully relaxed structure in the 
excited orbital state having Cih symmetry, and the en- 
ergy difference with the C^y configuration was 7.7 meV. 
The slightly larger energy difference for bigger supercell 
is expected, since more atoms around the defect can be 
relaxed. Note that the unstable modes have not been 
found in the recent calculations'^^, either due to different 
calculation method, or due to smaller size of the supercell 
used there. 

This symmetry breaking is not very surprising. The 
structure of the excited state corresponds to the typical 
case where the Jahn- Teller (JT) effect arises^"'^^: the sys- 
tem of C^y symmetry with two degenerate i?-type elec- 
tronic states and two degenerate e-type phonon modes 
{E®e case). In agreement with expectations, the energy 
difference between the parent high-symmetry Cst, config- 
uration and the symmetry-broken Cih structure comes 
mainly from the asymmetric relaxation of the atoms adja- 
cent to the vacancy. Therefore, although the magnitude 
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of the energy difference does not noticeably exceed the 
error of the DFT-PBE calculations, we expect that the 
Cih configuration is indeed the lowest-energy configura- 
tion at zero temperature within adiabatic approximation. 
This corresponds to the static JT distortion of the NV 
center's atomic structure, when the energy of the system 
has three potential energy minima, related to each other 
via rotation by the angle 27r/3 around the symmetry axis 
of the parent C^-u structure. At higher temperatures, ex- 
ceeding the C^y-Cih energy difference, (above ~ 100 K), 
all energy minima would be occupied equally, and an ef- 
fective symmetry would be observed in most (but 
not all) experimental situations'^. 

There is another possibility for the original sym- 
metry to be restored, the dynamical JT effect. Our calcu- 
lations have been performed within the adiabatic approx- 
imation, where the atoms are treated as classical point 
objects with well-defined coordinates. This approxima- 
tion may become invalid when two electron states are 
degenerate, as is the case for the C^y excited stated. 
The atomic motion then should be treated quantum- 
mechanically, and the vibrations may become strongly 
entangled with orbital degrees of freedom. When the JT 
energy barrier Ejt is much lower than the relevant vi- 
brational frequencies w, the relief of the potential energy 
surface is negligible in comparison with the zero-point 
atomic vibrations, and the total symmetry of the coupled 
electron-lattice system is C^y even at zero temperature. 
Recent experiments evidence in favor of the dynamical 
JT effect in NV centers^"^, while the static JT seems to 
contradict the existing datap^i*^. 

Our calculations seem to support the picture of the dy- 
namical JT effect, giving the barrier height Ejt of order 
of 10 meV, and the relevant vibrational frequencies be- 
ing an order of magnitude larger. This is the indication 
of the weak coupling regime, and is in agreement with 
the known experimental facts (such as e.g. small renor- 
malization of the g-factor, and the shape of the optical 
transition bands). However, detailed quantitative theory 
requires a special focused investigations, which would in- 
clude the lattice nonlinearity and the many-body effects 
associated with it. Such a study is to be performed in the 
future. Therefore, below we choose the Cih configuration 
to study the lattice vibrations in the excited state. For- 
tunately, the results below are not very sensitive to this 
choice. 



IV. LATTICE VIBRATIONS AND 
QUASI-LOCALIZED VIBRATIONAL MODES OF 
THE NV CENTER 

When a defect is introduced in a perfect lattice, the 
localized vibration modes may appear in the phonon 
spectrum^. If the energy of such a mode lies in the gap 
of the phonon spectrum, or above the highest phonon 
mode of the perfect lattice, then the vibration amplitude 
decays exponentially away from the defect (complete lo- 



calization). But when this mode lies inside the phonon 
spectrum of the perfect lattice, it can hybridize with the 
delocalized phonons, and give rise to a quasi-localized 
vibrational mode (qLVM) , which is not completely local- 
ized, but whose amplitude still concentrates in the im- 
mediate vicinity of the defect. The degree of the mode 
localization can be quantified by its inverse participation 
ratio (IPR), defined as 
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where N is the number of atoms in the system, and 
Uj is the normalized 3iV-dimensional vector of the 
atomic displacements of the corresponding vibrational 
eigenmode^. IPR characterizes a typical number of 
atoms involved in a given vibration, so that for com- 
pletely delocalized phonons P ^ 1/N is almost zero 
for macroscopic systems, while for completely localized 
modes P remains finite independently of the system size. 

The total density of states (DOS) of the lattice vibra- 
tions in the ground state of the NV center and the IPRs 
of the corresponding modes are shown in Fig.[3l The cor- 
responding graphs for the excited orbital state of the Cih 
symmetry are presented in Fig. |31 In both Figures, one 
can see a large number of modes with different degree 
of localization, including well-pronounced qLVMs with 
P of order of 0.03-0.05. Such modes are very strongly 
localized near the NV center, and involve noticeable vi- 
brations of only 4-6 atoms; the remaining contribution 
to the value of P is due to very small vibrations of more 
distant atoms. The list of modes with the largest IPR is 
given in Table U for the ground state and in Table [Tl] for 
the excited Cih state. Note that all modes are infrared- 
active. These results show that a large number of qLVMs, 
with different degree of localization, appears in the dia- 
mond lattice with the NV center. This is in contrast with 
the finding reported recently'^?-, where only six localized 
modes have been reported. This might be due a different 
computational method, or due to a smaller size of the 
supercell used in Ref. [S^. 

Upon the optical excitation, the qLVMs undergo sig- 
nificant changes, as clearly seen from comparison of 
Figs, [alb) and and Tables U and III The qLVM 

with the frequency of 41.66 meV (336 cm^^) which is 
clearly seen in the ground state, disappears in the ex- 
cited state, and a new qLVM at 31.62 meV (255 cm~^) 
appears in the excited state. Also, the group of qLVMs 
with frequencies slightly above 70 meV undergoes some 
changes: they become less localized in the excited state. 
Also, the group of the qLVMs at high frequencies, around 
140-160 meV, undergoes noticeable change. 

Note that most qLVMs are barely seen in the total 
DOS, Figs, ^a) and HJ^a), and only the IPR analysis 
allows us to discern them. To see the changes in the 
qLVMs in more detail, in Fig. [5]we present the projected 
DOS (PDOS) corresponding only to the vibrations of 
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TABLE I: Parameters of the quasi-localized vibrational modes with largest IPR in the Csu ground state: frequency u, IPR P, 
and symmetry. 

uj (meV) 41.66 71.21 72.16 74.79 97.95 138.63 156.64 156.70 

P 0.023 0.018 0.046 0.047 0.017 0.017 0.029 0.031 

Symmetry e e a a a e e a 



TABLE IL Parameters of the quasi-localized vibrational modes with largest IPR in the Cih excited state: frequency oj and 
IPR P. 



UJ (meV) 


31.62 


71.53 


73.42 


74.82 


97.70 


134.07 


156.85 


157.34 


162.69 


162.77 


163.35 


164.57 


P 


0.017 


0.060 


0.028 


0.034 


0.016 


0.017 


0.022 


0.028 


0.037 


0.017 


0.037 


0.082 




FIG. 3: (Color online). The total density of vibratioinal 
states (a) and the IPRs of the corresponding modes (b) for the 
ground state of the NV center. The graph (black solid line) 
in panel (a) was smeared, following the standard practice, by 
convolution of the calculated discrete set of frequencies with 
the Gaussian line of the width 0.05 meV. 



FIG. 4: (Color online). The total density of vibratioinal states 
(a) and the IPRs of the corresponding modes (b) for the ex- 
cited Cih state of the NV center. The curve (black solid line) 
in panel (a) was smeared, following the standard practice, by 
convolution of the calculated discrete set of frequencies with 
the Gaussian line of the width 0.05 meV. 



the nitrogen and the three carbons adjacent to the va- 
cancy. In the ground state, all carbons are equivalent, 
and have the same PDOS. In the Cih excited state, one 
carbon (denoted as Ci) is different from the two others 
(denoted as C2 and C3, and related to each other by re- 
flection in the mirror plane {110} passing through N, Ci, 
and the vacancy site), so that only the PDOS of the C2 
and C3 atoms are the same. Fig. [S] shows, in particular, 
that both the ground-state qLVM at 41.66 meV and the 
excited-state qLVM at 31.62 meV are associated almost 
completely with the carbon vibrations. The qLVM at 
41.66 meV is doubly degenerate (e type), so the contri- 
butions of diff'erent carbons are not well defined, while 
the excited-state qLVM at 31.62 meV is associated ex- 
clusively with Ci vibrations. 

Next, the group of qLVMs shghtly above 70 meV, 



which includes both nitrogen and carbon vibrations, 
shows moderate changes in the excited state. But 
the changes in the high-frequency qLVMs, around 140- 
160 meV, are very noticeable: the modes associated pri- 
marily with N vibrations shift towards lower energies, 
while the carbon vibrational modes acquire higher fre- 
quencies. These changes are likely to come from the elec- 
tron transfer from carbons to the nitrogen (from the D to 
the Ex orbital. Fig. [2]), which accompanies the optical ex- 
citation. Filling of the orbitals near carbon atoms makes 
the bonds stifFer, while the depletion of the nitrogen or- 
bitals makes them softer, with corresponding changes in 
the qLVM frequencies. 

Our calculations provide even more detailed picture of 
each relevant qLVM, identifying its type (bond stretch- 
ing, bond bending, etc.), and visualizing the motion of all 
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FIG. 5: (Color online), (a): Projected DOS (PDOS) of the 
vibrations associated with the carbon atoms adjacent to the 
vacancy. The ground-state PDOS of all carbons are the same 
(black solid line), while in the excited state the PDOS of the 
Ci carbon (dash-dotted green line) is different from the PDOS 
of the C2 and C3 atoms (dashed red line), (b): PDOS of the 
vibrations associated with the nitrogen atom. The ground- 
state PDOS is shown with black solid line, and the excited- 
state PDOS is shown with the dashed red line. All curves 
were smeared, following the standard practice, by convolution 
of the calculated discrete set of frequencies with the Gaussian 
line of the width 0.05 meV. 

atoms. In Figs. [Hand [71 we present the atomic vibration 
patterns of several important qLVMs for the ground and 
the excited state, respectively. In particular, the qLVMs 
with frequencies around 70 meV are shown: these modes 
are important for discussion of experiments given below. 
Detailed assessment of several other qLVMs is presented 
in the Appendix. Here we just mention that all the anal- 
ysis and the ideas presented above are consistent with 
detailed consideration of qLVMs. 

V. COMPARISON WITH THE 
EXPERIMENTAL RESULTS. 

The optical transition between the ground and the ex- 
cited state is accompanied by the changes in the atomic 
positions in the neighborhood of the defect site, i.e. in- 
volves absorption and emission of the lattice vibration 
quanta31i^. The zero-phonon line (ZPL) in the opti- 
cal spectrum, which corresponds to the vertical transi- 
tion (no absorption/emission of the vibrational quanta), 
is accompanied by sidebands corresponding to absorp- 
tion/emission of vibration quanta. For NV center, the 



FIG. 6: (Color online). qLVMs in the ground state of the NV 
center: atomic vibration patterns of several modes of differ- 
ent symmetry having largest IPR. The red arrows show the 
atomic vibrations, the arrow length is proportional to the vi- 
bration amplitude of a given atom, (a) and (b) — the modes 
of a symmetry, with energies of 72.16 meV and 74.79 meV, 
(c) — two degenerate modes of e symmetry with the energy 
of 41.66 meV appearing in the ground state but not in the 
excited state. Yellow spheres denote N atom, blue spheres 
denote the C atoms adjacent to the vacancy, green spheres de- 
note the C atoms adjacent to the N atom. Small grey speres 
are the other carbon atoms in the diamond lattice. 



ZPL has the energy 1.945 eV (640 nm), while the vibron- 
ically broadened sidebands extend well beyond 750 nm 
at room temperature. Since the optical transitions are 
tightly linked to the lattice vibrations, the studies of the 
optical emission and absorption provide clues about the 
relevant vibrational modes. In particular, in Ref. [sol the 
absorption and luminescence spectra have been studied 
in detail, and the DOS of the vibrations participating 
in the optical transition has been calculated. It has two 
broad peaks with maxima around 70 meV and 140 meV 
and width of order of 50 meV each, see Fig. [HI Later, 
more detailed examination of the PL spectr revealed 
the well-defined one-phonon sidebands, which correspond 
to a vibration mode of a symmetry with the energy 
of 70 meV. Moreover, the absorption band was differ- 
ent from the emission band, demonstrating at least two 
closely spaced peaks. 

Our calculations allow a direct comparison with these 
experimental results. First, let us focus of the DOS g{ijj) 
calculated from the experimental data in Ref. [s^. The 
qLVMs are expected to play primary role in optical tran- 
sition, in comparison with the delocalized phonons, and 
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FIG. 7: (Color online). qLVMs in the excited state of the 
NV center: atomic vibration patterns of several modes of dif- 
ferent symmetry having largest IPR. The red arrows show 
the atomic vibrations, the arrow length is proportional to the 
vibration amplitude of a given atom, (a) — the mode of fre- 
quency 31.62 meV appearing in the excited state, (b), (c), and 
(d) — the modes with frequencies 71.53 meV, 73.42 meV, and 
74.82 meV: these are the modes with largest IPR belonging 
to the group of qLVMs near 70 meV. Yellow spheres denote 
N atom, blue spheres denote the C atoms adjacent to the va- 
cancy, green spheres denote the C atoms adjacent to the N 
atom. Small grey spheres are the other carbon atoms in the 
diamond lattice. 



we expect that the experimentally extracted DOS is asso- 
ciated mostly with qLVMs, as is mostly governed by 
the vibration amplitude at the defect site^'*. In Fig.|Sl we 
compare the experimentally obtained giyi) curve (taken 
from Ref. Is^) with the total vibrational DOS, and with 
the projected DOS associated only with the vibrations 
of the nitrogen atom and the three carbons adjacent to 
the vacancy. For comparison with the experimental data, 
the discrete lines corresponding to individual modes have 
been broadened, following the standard practice, by con- 
volution with the Gaussian of the width 10 meV, cor- 
responding to the low resolution of the method used in 
Ref. [sO- The projected DOS practically coincides with 
the experimental data, thus confirming our expectation, 
while the total DOS demonstrates serious difference in 
the positions of the two peaks and in their relative am- 
phtude. Comparing Figs. |8] and [U we see that the two 
experimentally found peaks correspond to the two groups 
of qLVMs, one around 70 meV and the other around 140- 
160 meV. Note that the high-frequency group of qLVMs 
has not been reported before. These results clearly show 
the crucial influence of the qLVMs on the optical emission 
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FIG. 8: (Color online). The experimentally determined DOS 
of the vibrations affecting the optical transition (black solid 
line) are compared with the calculations results. The blue 
dash-dotted line denotes the calculated projected DOS of the 
vibrations associated with the four atoms (nitrogen and three 
carbons) comprising the core of the NV center. The red 
dashed line denotes the calculated total DOS of the lattice 
vibrations. An excellent quantitative agreement between the 
experiments and the calculated projected DOS is clearly visi- 
ble. Both calculated DOS were smeared by convolution with 
the Gaussian of the width 10 meV, corresponding to the low 
resolution of the experimentally determined DOS. 



and absorption. 

In the subsequent work'^^, the influence of vibrations 
on the emission/absorption has been analysed in much 
finer detail, at the level of individual vibronic peaks. It 
has been found that the vibronic spectra show a clear 
one-phonon sideband with the vibronic frequency of ap- 
proximately 70 meV. The experiments have been per- 
formed at low temperatures, so that the starting state 
for both absorption and emission is likely to contain no 
qLVM quanta, while the final state for absorption (emis- 
sion) will contain a single quantum of the qLVM cor- 
responding to the excited (ground) state^li^. We see 
immediately that the qLVM group with the frequencies 
near 70 meV, analyzed above and shown in Figs. IH] and 
[3 is exactly the right candidate for explanation of this 
feature. They also have correct a symmetry, see Tables U 
andHJ 

A puzzling feature found in Ref 'sT is that the absorp- 
tion sideband looks like a doublet of two individual lines, 
with the width of order of 20 meV each, and the peak- 
to-peak distance of only 3 meV (significant broadening is 
likely caused by the random stresses in the sample). Note 
that the multiplet structure of the one-phonon sidebands 
is fully consistent with our results: we have a group of at 
least three closely spaced qLVMs near 70 meV in both 
excited and ground state. The apparently single-peak 
sideband in emission and doublet in absorption, shown 
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in Fig. 6 of Ref. [Sll, are equally likely to be made of three 
peaks of different height, each inhomogeneously broad- 
ened by 20 meV. Therefore, we do not analyze this fea- 
ture in more detail. Possibly, future measurements on 
a single NV center, free of such a large inhomogeneous 
broadening, will convincingly resolve the fine structure of 
these sidebands. 

For a similar reason, we do not discuss the spectral sig- 
natures of the ground-state 41.66 meV and the excited- 
state 31.62 meV qLVMs. One might argue that the 
steeply rising edge of the one-phonon sideband in Fig. 6 
of Ref. [Sll does show a shoulder in the right position (near 
1.9 eV), but this claim would be way too speculative. 



VI. CONCLUSIONS 

We have performed detailed first-principles study of 
the localized distrortions and (quasi)localized vibration 
modes (qLVMs) of NV centers in diamond. Our calcula- 
tions show that the presence of NV center gives rise to a 
large number of qLVMs with different degree of localiza- 
tion. 

We found that in the high-symmetry excited state 
of the NV center some qLVM become unstable, and lead 
to lowering of the symmetry to Cih- Such a lowering 
is expected, since the excited state of the NV center is 
a typical example of Jahn- Teller unstable system. Our 
results seem to support the earlier suggestion'^'^ that the 
dynamical Jahn- Teller effect takes place. Our calcula- 
tions indicate the regime of the weak vibronic coupling, 
with the mode frequencies about an order of magnitude 
larger than the Jahn- Teller energy. While the quantita- 
tive first-principle theory of this effect for the NV center 
is lacking, it constitutes an exciting topic for future work. 

Next, we have analyzed and described some interest- 
ing qLVMs in the ground and excited orbital states of NV 
center, and discussed the changes in the qLVM density 
of states occuring upon optical excitation. We have per- 
formed a comparison of our results with the experimental 
data reported in Refs. lsoUsil We have shown that the cal- 
culated DOS of the qLVMs is in very good quantitative 
agreement with the results of Refs. [30|, demonstrating 
two peaks near 70 meV and 140-160 meV. These peaks 
are determined by the corresponding groups of highly lo- 
calized modes. Therefore, as expected, qLVMs crucially 
impact the optical properties of the NV centers in dia- 
mond. 

We have also demonstrated that the group of qLVMs 
with frequencies near 70 meV is likely to explain the fea- 
tures of the one-phonon sideband observed in Ref. IsTI . 
possibly including their multiplet structure. However, 
detailed resolution of this issue requires further exper- 
imental and theoretical work, such as high-resolution 
measurements of the absorption/emission spectrum on a 
single NV center, and highly precise first-principles cal- 
culations of the electron-phonon coupling. 
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Appendix A: Atomic displacement patterns of 
several ground-state and excited-state qLVMs 

Below, we present the patterns of atomic displacements 
for several ground-state and excited-state qLVMs having 
larges IPR. 

In the ground state, the lowest-energy localized a 
mode at 72.16 meV is shown in Fig. [Sl^a). It is strongly 
localized on the nitrogen: the N atom moves in the di- 
rection perpendicular to the plane defined by the three 
C atoms around the vacancy. The two a modes shown in 
Fig. [HJb) and (c), with the frequencies 74.79 meV and 
97.95 meV are the breathing modes of the four first- 
neighbor atoms of the vacancy: all four atoms move to (or 
away from) the vacancy. The displacements of the four 
first neighbor atoms in the 97.95 meV mode are smaller 
than those in 73.79 meV mode indicating less localization 
of the breathing mode at 97.95 meV. The other a mode 
at 156.7 meV, shown in Fig. EJd), is a twisting mode 
around the NV center. The first- and second-neighbor 
carbon atoms from the vacancy rotate in different direc- 
tions. 

The eigenvectors of the four pairs of doubly degenerate 
e modes, at the frequencies of 41.66 meV, 71.21 meV, 
138.63 meV, and 156.64 meV, are shown in Fig.|ni[e)-(h). 
The two lower-energy localized e modes. Fig. [9je) and 
(f), are mostly C-N-C angle bending modes. The third 
one, shown in Fig. Hfg), with a smaller IPR, is the N- 
C stretching mode at the NV center. The qLVM shown 
in Fig. IHJh) is a C-C stretching mode: carbon atoms 
around the vacancy move in the out-of-plane direction, 
and at the same time, their neighboring carbon atoms 
move in the opposite direction. 

The displacement patterns of the qLVMs in the Cih 
excited state are shown in Fig. I10[ and are quite differ- 
ent from the ground-state qLVMs. There is no doubly- 
degenerate mode in the structure with C'lh symmetry. 
Among the twelve localized modes with largest IPR, only 
one, with frequency 71.53 meV, shown in Fig. [TOlb). 
is similar to the ground-state mode at 72.16 meV. 
Six qLVMs at 134.07 meV, 157.34 meV, 162.69 meV, 
162.77 meV, 163.35 meV, and 164.57 meV, shown in 
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Fig. fTOTf.h.i.i.k.l') respectively, are the N-C and C-C 
stretching modes. There are three locahzed C-N-C angle 
bending modes at 31.62 meV, Fig.HU^a), at 74.82 meV, 
Fig.lHd), and at 156.85 meV Fig.HUIg). One breathing 
mode is seen at 97.70 meV, Fig.fTOle). which has smaller 
IPR. Finally, the mode at 73.42 meV has a C atom mov- 
ing opposite to the vacancy, as shown in Fig. lTUT c). In all 
the excited-state qLVMs discussed above, the Ci carbon 
atom moves independently, while the other two equiva- 
lent carbons move collectively to maintain the Cih sym- 
metry. The carbon atoms adjacent to the displaced car- 
bons move in the same plane. This plane is formed by 
the three first-neighbor carbons in such a way that the 
rotation and translation of the whole system is absent. 
Due to the low Cih symmetry, all these modes are IR 
active. 



FIG. 9: (Color online). qLVMs in the ground state of the NV 
center: atomic vibration patterns of several modes of different 
symmetry having largest IPR. Yellow spheres denote N atom, 
blue spheres denote the C atoms adjacent to the vacancy, 
green spheres denote the C atoms adjacent to the N atom. 
Small grey speres are the other carbon atoms in the diamond 
lattice. 
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